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The inclusive structure function ^2(0;, Q^) has been determined with extraordinary 
accuracy down to very smaU x by recent experiments at the HERA coUider [|l|. For 
Q^ ^ IGeV the x and Q^ dependence of the data is in complete agreement with that 
predicted by the next-to-leading order (NLO) Altarelli-Parisi evolution equations. At small 
X and large Q^, the evolution equations are dominated by the small x singularities of the 
Altarelli-Parisi splitting functions ||^, and retaining only the singularities in the LO and 
NLO splitting functions yields an excellent approximation to the full solution in the HERA 
region |^. 

As we go to higher orders in ctg the splitting functions become more and more singular, 
and these higher order singularities might be expected to become dominant at small enough 
X. It thus appears reasonable to try to improve the description of small x evolution by 
supplementing the usual leading-order splitting functions with contributions which sum all 
leading logs of x (LLx) to all orders in a^, i.e. all terms of the form (a^ logx)"^. Likewise, 
the NLO splitting functions can be supplemented by a summation of next-to-leading log 
X contributions {i.e. all terms of the form as(as logx)"', and so on (the 'double-leading 
expansion' |^). 

However, as is by now well known [^, such attempts are unsuccessful: essentially, the 
data in the HERA region are so accurately described by plain NLO evolution (and the 
small X approximation to it) that any further "improvement" would spoil this agreement 
unless its effects were extremely small. Furthermore, the recent determination [§|,0J^] of the 
subleading corrections to LLx evolution has shown that NLLx contributions are extremely 
large, and in fact grow faster as x — > than their LLx counterparts [^,1C]. Hence, contrary 



to naive expectations, the double leading expansion does not appear to be stable at small 
X. A deeper understanding of the all-order behaviour of splitting functions at small x is 
required. 

In this letter we will determine the small x behaviour of the Altarelli-Parisi splitting 
functions order by order in LLx, NLLx,. . . . All our discussion will be based on a formal per- 
turbative treatment of the small x expansion: at NLLx we retain only the contributions to 
small X evolution calculated in [^ , systematically discarding any terms which are formally 
NNLLx. This approach allows us to isolate the reason for the asymptotic breakdown of the 
small X expansion. This turns out to be unrelated to various problems discussed elsewhere, 
such as the unphysical behaviour of the solutions found in [|ll|,|l2| and the running coupling 



resummation effects discussed in [p!3| , p^ . We find that the instability observed at NLLx is 
not some peculiar feature of the calculation of ref. , but is completely generic, probably 
persisting to all orders. Formally subleading contributions to the splitting function are 
not suppressed by powers of a^, but rather grow faster and faster at small x. We derive 
conditions on the small x evolution kernel which are necessary for stable evolution, and 
show that they may be satisfied by a suitable reorganization of the perturbative expansion. 
Such a resummation introduces an a priori undetermined parameter which describes the 
resummed all-order small-x behaviour of the structure functions. 

The small x behaviour of splitting functions P{x) is most easily studied by considering 
their Mellin transforms, the anomalous dimensions 7(A^) = Ldxx^ P{x). The leading 
small X behaviour is then found by expanding 7(A^) about its rightmost singularity in the 
complex A^ plane, which in the singlet sector is located at A^ = (and at A^ = — 1 in 

1 



the nonsinglet sector, which is thus suppressed by a power of x and hence asymptoticaUy 
neghgible). The anomalous dimensions in the singlet sector are given by a two- by- two 
matrix; however only one of the two eigenvalues of the matrix is singular at A^ = at 
leading order (and to all orders in appropriate factorization schemes). It is thus sufficient 
to concentrate on this leading eigenvalue and its associated eigenvector Gn{Q'^), the Mellin 
transform of the distribution function G{x,Q'^), which satisfies the evolution equation 

j^G^iQ')=^iN;a)GNiQ'), (1) 

where t = ln((5^/A^) and for future convenience we write a{t) = as{t)/2n. 

The general structure of the anomalous dimension in the small x expansion is 

7(iV;a) = 7o(a/iV) +a7i(a/iV) + ... (2) 

-'(^)-E4"'(8iV.ln2^)". (3) 

n=l 

where the coefficents A)^ have been normalized such that the radius of convergence of the 
series is one. The associated splitting functions Pk are immediately obtained by inverse 
Mellin transformation of 7^ : 

P{x;a) = -[Po{E)+aP,{E) + ...] (4) 

X 









n=l 

where ^ = In - and we define 

E = (2CA)(41n2)^e = 8NJn2aC (6) 

Subsequent terms 7^, Pk in the expansions ®,(^ of the anomalous dimension and 
splitting function sum the leading, next-to-leading,. . . logarithms of -. They can therefore 



be determined [ p!5| , p!6[] from knowledge of the respective leading, next-to-leading,. . . QCD 
high-energy asymptotics, as given by leading log x evolution, which is in turn controlled 
by an equation of the form 

—GM{x)=ax{M;a)GM{x), (7) 

where the Mellin variable M is defined by 

Gm(x) = ^°°^(|J) G(x,Q2) (8) 

and the anomalous dimension x(M; a) admits a perturbative expansion 

x{M-a) = xo{M) + axiiM) + ... . (9) 



The leading-order term XoiM) = 2Nc[2il;{l) - ^{M) - ^{l - M)] is well known fl^, 
while the next-to-leading term xi has been determined only recently [^.q Note that since 
structure functions scale in the Q^ -^ oo limit and drop linearly with Q^ as Q^ — > the 
leading- twist physical region corresponds to < M < 1. 

The /c-th order small x anomalous dimensions 7A; (;§^) eq. (^ can be determined from 
Xo, • • • 5 Xfc eq. (^) by matching the solutions to the respective evolution equations. Assume 
for the moment that the coupling is fixed. Solving eq. (|^) by Mellin transform with respect 
to A^, we find 

where G^{M) is the boundary condition at ^ = 0. In order to compare to eq. (|1|), invert 
the M- Mellin at large Q^: 

where Mp{N; a) is the position of the rightmost pole of Gnm in the M-plane in the physical 
region < M < 1. It follows that the solutions to the evolution equations (|l|) and (]^) 



coincide only if 7 = M„, i.e. if their anomalous dimensions satisfy the 'duality' relation WB 



N 
x[7(A;a);a] = -. (12) 

a 



Expansion of eq. (12) in powers of a keeping a/N fixed gives a set of relations which 



determine 7/c perturbatively order by order in terms of the perturbative expansion of x'- 

N 
Xoho{a/N)] = - (13) 

7i(a/A)^-2aM^ (14) 

Xono{a/N)\ 



X2(xo')' - XiXi^Xo' + |(xi)^Xo" 
(Xo' 



72(a/A) = ;-— 3 , (15) 



and so forth, where the prime indicates differentiation with respect to M, and in the last 
equation all x- functions have argument 70 (a/A). Eq. ([T^ ) should be viewed as an implicit 
equation for 70, with xo{M) evaluated in the physical region < M < 1. 

The derivation presented so far holds at fixed coupling a^. If the coupling runs with Q^ 
in both equations, it is sufficient to include the running up to order ct^ when computing the 
anomalous dimensions to k-th order in the LLx expansion. The contributions to the kernel 



Note that here we have adopted a different set of normalization conventions to those used 
in ref. [P, where the right hand side of eq. (^ is written as 2Nc{x + ^l""' ^)'- o^^ Xo is the same 
as their 2NcXj while our xi is the same as their N^S- 



on the right hand side of eq.(|^ are then found by replacing a^ by differential operators: 
for example, to NLLx it is sufficient to let 

a^a(l + 6oa^j, (16) 

where bo = ^{^Nc — ^Uf). Solving (^ as before, taking care to treat all NLLx terms 
(including the derivative term) as perturbations, now gives 
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Inverting the M-Mellin, and again comparing with the solution of eq. (|ip with the running 
coupling expanded at NLLx, we see that the solutions match if in place of ([T^ ) we now 
have 

f /An 1 / / N , u Xo(7o)Xo(7o)\ .^gv 

V-Xo(7o), (19) 



Xi(7o) 1 d 



Xo'(7o) adt 



In 



since to leading order da/dt = —hoo?, while axo(7o) = A^ so d'-yo/dlna = — Xo/Xo- The 
expression (^U[) was obtained in |T^ by solving eq.(|^ with running coupling exactly, and 
then inverting the Mellin by a saddle point argument consistently at NLLx. 

When the coupling runs there is however a further ambiguity related to the choice of 
factorization scheme. Under a LLx change in the normalization of Gm {0 5 ^- ^- ^m (0 ~^ 
u{M)Gm{0^ tti6 NLLx kernel changes according to 

XI (M) ^ XI (M) - 60X0 (M) -^ In u{M) . (20) 

From eq. (p!9|) the anomalous dimension changes as 

7^7 + — lnw(7o). (21) 

This is the same as the NLLx shift in the anomalous dimension induced by the LLx scheme 
change GN{t) -^ u{a{t)/N)GN{t), where u{a/N) = 1 + ^f ^^(a/iV)'" 0. Note that a 



NLLx scheme change only affects the anomalous dimension at NNLLx: the mismatch is 
due to the coupling running with Q^ while the logarithms are ordered in x. It follows 
that knowledge of the leading order coefficient function [^0[ is sufficient for a consistent 
calculation at NLLx. 

In order to determine the NLLx anomalous dimension 71 in MS factorizationH from the 
Fadin-Lipatov kernel xf^ (eqns.(14)& (22) of ref.0) two further adjustments are required. 



Note that the large eigenvalue of evolution at small x is invariant under the usual (NLLQ) 
scheme changes: it fohows that for our purposes the MS scheme is the same as the DIS scheme, 
and indeed we will use these two notations interchangeably. 
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First, the kernel xi for the smaU x evolution of the distribution Gm{0 — ^~^9M{i) 
is related to that for the unintegrated distribution QmH) employed in [^ by xi(M) = 
xf^(M) +6oXo(-^)/-^- This may be thought of as a LLx scheme change with u{M) = M. 
Second, the correct expression for 71 with MS factorization requires a further scheme 
change with u{M) = R{M), where i2(7o(a/A^))) = RNi^s) is calculated in ref.|2g].i The 
MS NLLx anomalous dimension is thus 



7.(a/iV)^-^ + i|lni.(7o)7oV-X^(7o) 



(22) 



Xo'(7o) 



xr (70) + 6oiVc((2^'(l) - ^'(70) - ^'(1 - 70)) + 4^X0(70)') , (23) 



where in the second line we used eq.(B.18) of [20] for dlnR/d'fo. 

It is clear from eqns.(^) & ( ^Bj ) that at NLLx the effect of the running of a^, eq.([TPD, 
can be entirely absorbed into a shift of the NLLx anomalous dimension 71 equivalent to a 
choice of factorization scheme. It turns out that this property persits to higher orders in 
the small x expansion. We may thus view running coupling effects as a contribution to x'- 
henceforth we assume that x incorporates the running coupling effects at the appropriate 
order in the small x expansion, so that the anomalous dimension 7 is given by the duality 
relation (|T2l). 



To see this compare the result for 7^^ in MS obtained in ||2^ with that derived in |21]: their 
ratio must be u. Alternatively, the explicit derivation given in Appendix B of pC[| of the anomalous 
dimension from the dimensionally regularized kernel may be extended to NLLx: this gives the 



same result [22|. 




Figure 2: The LLx splitting function Po(S) (solid, positive) and the NLLx splitting 
function ^Pi{E) (negative), computed numerically (using eqns. (p|),(p^,(19),([2^) and 
Xi^ from ref. [^), and renormalised by a factor e^", in various factorization schemes: 
DIS (or MS, see footnote 2) (solid), Qo-DIS (dotted), SDIS (23| (dashed) and GDIS 
MM (dot-dashed). 



As is well known, the BFKL kernel xo(-^) is symmetric about M = i, where it 
has a minimum, with Xo(|) = 8A^cln2. It follows that 70, viewed as a series of powers 
of 8Nc\n2-^ (eq.(^), has unit radius of convergence. All higher order 7/c will have the 
same radius of convergence in any factorization scheme in which all Xk{M) are free of 
singularities for < M < i. The fact that the series for the anomalous dimension has 
finite radius of convergence is important because it shows that the corresponding splitting 
function has infinite radius of convergence as a power series in S, and can thus be used 
down to arbitrarily small x, provided a large enough number of terms is included Q. In 
the asymptotic region where Aj^ /^k ^ ^ ^^^ number of terms which must be included 
is of order kc{x^ Q^)^ where kc is the solution of the implicit equation 'E^" /k^^ ~ 1, whence 
kc ~ 2.7E] for S ^ 1. It is apparent that already in the HERA region, where S can be as 
large as ten (see fig. 1) a large (though finite) number of terms should be included.a 

Since each contribution Pk to the small x expansion (^) of the splitting function sums 
up leading logs of x to all orders, each subsequent order in the expansion appears to be 
of order a^ compared to the previous order. It was thus natural to conjecture ||T3| , pU| ,^ 
that evolution at small x can be accurately described by truncating the expansion to finite 
order, provided only that 0;^ is small enough {i.e. Q^ is large enough). As discussed in the 



It is interesting to contrast this with the familiar case of large x (Sudakov) resummation, 
where the series of contributions to be resummed in the anomalous dimension is divergent and 
can only be treated by Borel resummation: for meaningful results an infinite number of terms 
must be included whenever as ln(l — x) ~ 1. 
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Figure 3: The ratio ^Pi{E)/Po{E) of NLLx and LLx splitting functions, computed as 
in fig. 2, in tlie same factorization schemes . 



introduction, however, this conjecture does not seem in agreement with phenomenology 
nor with the numerical comparison of the NLLx splitting function Pi to the LLx Pq MKO 



as may be seen from Fig. 2, Pi is large and negative, and furthermore the size of the ratio 
of Pi to Pq (see fig. 3) increases rapidly at small x, i.e. as S increases. It is apparent from 
the plot that the ratio ^Pi/Pq is much greater than l/as{Q'^) ^ 10 throughout most of 
the range of S relevant for HERA and the LHC. This in itself explains the failure of the 
phenomenology based on a leading-order truncation of the expansion at HERA: the small 
X expansion breaks down for any reasonable value of ctg [Q . 

To understand this result we will now consider analytically the asymptotic behaviour 
of Pk at large values of S. The asymptotic behaviour of the splitting functions may 
be determined through their definition as the inverse Mellin transform of the anomalous 
dimension 7 [^: at LLx 



Pn(E) 



dN 
2nia 



.^N 



7o{a/N) 



, ZTTZ 



(24) 



where we have used eq. (^3|) to change integration variable dN — aXo{7o)d'lo- The asymp- 
totic expansion of Po(S) as ^ — ;> 00 can then be determined by the saddle point method: 
a straightforward computation, remembering that the only real minimum of Xo(7o) in the 
range < 70 < 1 is at 70 = | leads to 



^o(S) 



X0C2) 



1/2 



Xo(i) = 

— p 

^3/2 ^ 



1 + 



xo(i)xr(^)i^^,^ 



8(X[,'(^)) 



1 



(25) 



(n)fV 



where S is given by eq. (P), while the n-th derivatives Xo (l) ~ 'iNcn\{2'^~^^ — 1)C(^ + 1 
when n is even (all odd derivatives vanish). 



The asymptotic behaviour of the next-to-leading correction Pi can be determined in 
a similar way: 



Pii^) 



°° dN 



.^N 



Xiboja/N)] 
x'obo{a/N)] 

1/2 /I 



drfo_ 
2tii 



>€axo(7o) 



Xi(7o) 



-2''E 



Sl/2 



i + o - 



(26) 



where in the last step we have assumed xi(M) to be regular at M = ^ 



This require- 



ment can always be achieved by choice of factorization scheme: in particular Xi^{\) 
-71.64A^2 _ o.52A^cn/ - IQ.luf/Nc is finite [§1 and in the MS factorization scheme xi{ 
is effectively (see eq.(|2|)) xf^{\) + hoN^[{A\iv2f - ^]. 

It follows that asymptotically at large S in the MS scheme 

^i(S) Xi{\)^ 



Po(S) -^^ xo(i) 



S + 0(1) 



(27) 



the next-to-leading correction, despite being suppressed by a factor of ctg, rises linearly 
with S and hence with Inl/x at small x, and becomes eventually dominant. In terms of 
the anomalous dimension 7, this means that the ratio A)^ /A^' of the NLLx coefficients 
to the LLx coefficients in the expansions (|^) rises linearly with n. It is clear that the origin 
of the rise is the simple pole of 71 eq. ([T^), viewed as a function of 70, at 70 = |. Because 
the denominator of eq. (|I4D vanishes linearly at 70, this simple pole is present whenever 
Xi{\) has a finite nonzero value. 

In factorization schemes where xi(M) diverges at M = 1/2 the singularity in 71 will be 
stronger, and so Pi will rise more rapidly at large ^. Such factorization schemes have been 
considered in the literature as being possibly more appropriate at small x: in particular, 
the Qo-schemes [0, the SDIS scheme Q and the "physical" scheme (GDIS) 0|li, all of 
which have the advantage of reducing the size of the leading perturbative corrections in the 
quark sector. This reduction is accomplished at the expense of introducing a singularity 
in xi(Af) at M = i: for instance in the Qo^scheme (i.e. eq.(|T^)) xi(M) has a simple pole 
at M = ^ due to the vanishing of Xo ( ^ ) • 

It is easy to see that if the NLO anomalous dimension has a simple pole at the location 
of the LO saddle, this dominates the asymptotic behaviour of integral (pO) . In this case. 



on top of the contribution of eq. (|2q ), there is a further contribution from the residue of 
the pole: 

Pr^-(S) ~ -iRes[xi]e=-f..., (28) 

=L — »-00 

where Res[xi] is the residue of the simple pole of xi(-^) at M = i: it is equal to |&oXo(^) 
in all of the singular schemes mentioned above. The ratio Pi/Pq now rises asH S"^'^: 

1/2 



A(S) 
i^o(S) 



Xi 



'H\) 



Xo(i) 



+ /c S - 60 



8xo(^) 



73/2 



1 



xo(i)xr(^) i 

8(X[,'(^))2 S, 



+ 0(1), (29) 



Note that the results given in both [P] 
A similar result was found in ref. |13]. 



and [p5| are numerically incorrect. 




Figure 4: As fig. 3, but with the asymptotic behaviours ( |27] ) and 
lytically subtracted from the previous results computed numerically, 
in scale on the vertical axis. 



) calculated ana- 
Note the change 



where /c = 26o, -26o, in the Qo-DIS g, SDIS Q and GDIS 12|,|T9| schemes respectively. 
In all such schemes, the small x expansion appears thus to be particularly badly behaved. 

In Fig. 4 the splitting function ratio Pi/Pq is plotted again but now with the asymp- 
totic results ( |27| ) and (^9|) subtracted. It is clear that the asymptotic behaviour sets in 
surprisingly quickly: already at S ^ 3 the subtracted ratio becomes constant. It follows 
that what is left of Pi after the subtraction is no longer unnaturally large: the asymptotic 
growth of the ratios ( pT]) and (^9[) is entirely responsible for the breakdown of small x per- 
turbation theory at NLLx. This could have been anticipated even before the calculation 
of 0: the NLLx correction to the LLx splitting function will inevitably become large 
asymptotically unless the NLLx correction to x(^) vanishes. 

Pursuing the argument to higher orders, it is apparent that the anomalous dimensions 
7i have higher order poles as the order of the expansion increases: 72 eq. (^3|) has a 
triple pole, and in general 7^ has a (2n — l)-th order pole. Consequently, the associated 
splitting functions will display stronger and stronger rises with ^. It follows that the 
small X expansion eqns. (0),(§) inevitably breaks down at small x: as S grows, the higher 
orders of the expansion become more and more important. This means that the leading 
contributions at small x have not been properly resummed. 

The origin of this failure can be simply understood by recalling that at fixed coupling 
the leading asymptotic small x behaviour of the solution to the ^-evolution eq. (|^) is given 
by x'^^^^^''' where Mg is the position of the saddle point in M (so at leading order x = Xo 
and Ms = ^)- When solving the small x Altarelli-Parisi equation eq. (P, this growth at 
small X, rather than being generated by solving an evolution equation in ^ eq. (|^) (with ^- 
independent anomalous dimensions), is included in the splitting functions. But expanding 
the corrections to the LLx asymptotic behaviour in powers of ag 



X 



a(xo(M,)+axi(M,) + ...) 



X 



aXo(2) 



Xo(2) 



(30) 



it is clear that the LLx asymptotic behaviour can be modified by subleading terms only 
if these rise with S. For instance, the NLLx correction can only be generated if Pi/Pq 
rises linearly with S, with slope Xi(|)/xo(|)5 as indeed we found above in eq. ([27|). At 
higher orders in a^ the LLx behaviour receives corrections proportional to higher powers 
of S, and correspondingly the higher order P^ have higher order poles. However these 
corrections are no longer given by a trivial exponentiation of the NLLx result. 

The bad behaviour of the small x expansion is thus due to the fact that the asymptotic 
behaviour of the solution to the evolution equations at large S does not coincide with the 
LLx prediction. The subsequent mismatch in the order of subleading corrections makes a 
nonsense of the perturbative expansion (^ and (|5|) . This can only be corrected by suitably 
reorganizing the expansion (^ in order to properly resum the large corrections. 

Since a change in the factorization scheme mixes different orders in the perturbative 
expansion, different scheme choices may be thought of as resummations. Furthermore 
just as there are scheme choices which make the perturbative expansion less stable, so 
there are choices which can improve it. More precisely, we can resum the large corrections 



by choosing the scheme in such a way that the anomalous dimensions eqns. (|T4[) , (p!5| ) , . . . 
are all regular at 70 = | order by order. This is always possible because if all 7^ with 
^ < ^0 ~ 1 are regular at 70 = ^, then 7^^ has a simple pole at 70 = | (assuming that 
all Xi (-^) are regular there) , which can be removed by choosing a scheme which subtracts 
a constant from Xi{M) equal to the numerators of eqns. (|np,([T^) and their higher order 
generalizations. Necessary conditions for a satisfactory perturbative scheme are thus that 
in the new scheme 

Xi(i) = 0, X2(i) = i(x;(i))Vxo(^), ■■•• (31) 

Clearly these conditions are not very restrictive. In fact, since there is only one 
condition at each perturbative order, they can be imposed simply by a choice of renor- 
malization scale, i.e. by the replacement of as(Q^) with ^^((/cQ)^), where k may itself be 
expanded as a series in a^. Then for example at NLLx A]^' -^ A]^' + A^'ubo log A;^, and 
the linear rise of A]^ /Aq in MS factorization schemes may be eliminated by choosing 
60 log /c = — ^Xi(^)/Xo(^)5 which gives k ~ 300. Choosing such a large scale does indeed 
lead to stable perturbative behaviour (see fig. 5). However it is also clearly a fine tuning: 
varying the scale by a factor of two either side leads to huge variations in the relative size of 
the perturbative correction. Other scale choices, such as BLM, designed to reduce the size 
of subleading corrections have been considered in ref. []5^ .0 At NLLx it is even possible to 
fine tune the choice of scheme such that xi(-^) = 0, by choosing the gluon normalization 
factor u{M) as a solution to the first order differential equation (Inu)' = xi{M) /boXoiM) 
(cf. (pUD): then 'yi{as/N) = and all NLLx corrections have been removed from the large 
eigenvalue of evolution. 



It is also possible to remove the singularity (28) in singular schemes such as Qo-DIS by tuning 
the scale. However now the choice of scale depends on x: A; ~ ki exp(/c2\/C + fca/y^), where each 
of fci, /c2 and fca all require fine tuning if the perturbative expansion is to be stable. Such a scale 
might also be justified through a BLM prescription 

10 



The physical meaning of these scheme choices is that, after the scheme change, the 
large corrections to the LLx asymptotic small x behaviour are absorbed into the x- 
dependent initial condition to the perturbative evolution in Q^. However because these 
large corrections have a large scheme dependence, it seems pointless to consider them at 
any finite order: all contributions to the asymptotic behaviour should be resummed to all 
orders, and then included in the LLx anomalous dimension. To this purpose it is sufficient 
to subtract Q'^-independent contributions in eq. (H): 

X(M; a) = c{a) + xo(M) + axi(M) + . . . , (32) 

where Xi(M) = Xi(M) — q, i = 0,1,2,... and formally c{a) = cq + J2'^=i '^nd^ ■ The 
subtractions ci, C2, . . . are then fixed by the criteria (^l]): we need 



Cl 



Xi(i), C2 = X2(i)-i(x;(i))Vxo(^), ■•■• (33) 



Since the resummed splitting function will be independent of cq, there is in principle no 
need to fix its value; however it is convenient to choose cq = xo(l/2) as with this choice 
the parameter 

A = ac{a) (34) 

has a direct physical meaning. 

We can now use the expansion (|3^ ) of x ^^ 6q.(^2]) to determine 7^ order by order, 
treating c{a) as leading order. At LLx the 'resummed' anomalous dimension 70 (A^; a) is 
then the solution of 

cia)+xo[lo]=N/a. (35) 

This implies that ^o{N; a) = 'yo[a/{N — {X — aco))] and consequently that the LLX splitting 
function 

Po(x;a) = Po(S)e(^-«^o)«. (36) 

In particular the asymptotic behaviour (|2^) becomes 






1 + - 



(37) 



The parameter A thus determines the nature of the asymptotic small x behaviour. The 
expansion in powers of a^ is now well-behaved: due to the conditions (^) all higher order 
Pi{x) behave in the same way as Pq at large ^, and thus all the corrections to Pq ^^^ down 
by powers of as uniformly in x. The ratio Pi/Pq is shown in fig. 5: it is indeed uniformly 
bounded and not too unreasonably large. 



In fact it is identical to the subtracted ratio in DIS plotted in fig. 4: this is because 

by an integration by parts and change of variables. 
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Figure 5: As fig. 3, but after various resummations: the subtraction ( |32| 
scale (dashed), and scales a factor of two either side of it (dotted). 



the fine tuned 



The reorganization of the perturbative expansion (|3^) can be viewed as an effective 
resummation of the higher orders of the expansion. Since c{a) — cq is of 0(a), eq. ([3^ ) 
imphes that Pq and Pq differ by a series of formally subleading contributions. However, 
the parameter A (^4|) which summarizes the asymptotic behaviour at small x is treated as 
Os-independent, and thus included in the leading order eq.(|35|). This 'order transmutation' 
effectively resums into the LLx anomalous dimension the all-order behaviour as given by 
A. However, to determine the value of A it may be necessary to use arguments which go 
beyond mass factorised perturbation theory. In particular, the unitarity constraint would 
suggest that A < 0: if A were positive the resulting powerlike growth in the splitting 
function at small x would drive a corresponding rise in the cross-section, which would 
ultimately violate the Proissart bound. 

The removal of the unbounded growth of formally subleading corrections at small x, 
achieved by the resummation described above, is a necessary prerequisite for a consistent 
small X resummation of Altarelli-Parisi evolution. Although our resummation does not 
resolve the instability in the small x evolution equation discussed in [ pTj - p^ , it does show 
that this is a separate issue. Indeed, the instability is clearly related to the shape of x as 
M — !► and M ^ 1, whereas the resummation criteria (|3T| ) refer to M = ^. At small M 
(and thus large Q^), the relevant approximation is to use the conventional Altarelli-Parisi 
equation: from this it may readily be inferred (using a duality argument [16| inverse to 
that used to obtain (|12|)) that the resummed kernel must always be finite and positive at 
M = 0, which is probably sufficient to cure the instability. Possibly related attempts to 
deal with these instabilities have been presented in ref. [^ . 

To conclude, we have shown that the poor behaviour of the small x expansion which is 

m 



manifested 
determination 



in the NLLx splitting functions computed from the recent Fadin-Lipatov 
of the next-to-leading high energy QCD asymptotics can be traced to 



the fact that the formally NLLx corrections to the LLx contributions to the splitting 
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functions are not truly subleading at small x. We have shown that this problem persists 
to all orders, and is related to the fact that the leading small x behaviour is not given by 
the leading order term of the small x expansion, but rather must come from an all-order 
resummation. We have demonstrated that a reorganization of the perturbative expansion 
is necessary, and given criteria eq. (^lD which must be met if such a resummation is to 
be successful. We further constructed a resummation which meets these criteria, but 



depends on a new parameter A eg. (|3^) which controls the asymptotic growth at small x. 
A complete resummation of Altarelli-Parisi splitting functions at small x might now be 
achieved through careful matching in the high Q^ region. 
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